Cooling toolbox for atoms in optical lattices 
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We propose and analyze several schemes for cooling bosonic and fermionic atoms in an optical 
lattice potential close to the ground state of the no-tunnelling regime. Some of the protocols rely 
on the concept of algorithmic cooling, which combines occupation number filtering with ideas from 
ensemble quantum computation. We also design algorithms that create an ensemble of defect-free 
quantum registers. We study the efficiency of our protocols for realistic temperatures and in the 
presence of a harmonic confinement. We also propose an incoherent physical implementation of 
filtering which can be operated in a continuous way. 

PACS numbers: 03.75.Lm, 03.75.Hh, 03.75.Gg 



INTRODUCTION 



Ultracold atoms stored in optical lattices can be con- 
trolled and manipulated with a very high degree of pre- 
cision and flexibility. This places them among the most 
promising candidates for implementing quantum compu- 
tations [E 13, S S 13 and quantum simulations of certain 
classes of quantum many-body systems 0, 0, S 0, EU 
ITU IliL ll3L IhH Il5j . Quantum simulation would allow 
us to understand physical properties of certain materials 
at low temperatures that so far have eluded a theoreti- 
cal description or numerical simulation. However, both 
quantum simulation and quantum computation with this 
system face a crucial problem: the temperature in cur- 
rent experiments is too high. In this paper we propose 
and analyze several methods to decrease the temperature 
of atoms in optical lattices and thus to reach the interest- 
ing regimes in quantum simulations, as well as to prepare 
defect-free registers for quantum computation. 

So far, several experimental groups have been able to 
load bosonic or fermionic atoms in optical lattices and 
reach the strong interaction regime \WL IT^ . ITsl \l$L I20L 
I2TL 12^ . I2H |24| . The analysis of experiments in the Tonks 
gas regime indicates a temperature of the order of the 
width of the lowest Bloch band an d for a Mott In- 
sulator (MI) a temperature of the order of the on-site in- 
teraction energy has been reported 0, |2!| . For fermions 
one observes temperatures of the order of the Fermi en- 
ergy j2|| HU- Those temperatures severely restrict 
the physical phenomena that can be observed and the 
quantum information tasks that can be carried out with 
these lattices. 

One may think of several ways of cooling atoms in 
optical lattices. Since the process of loading atoms into 
the lattice may lead to additional heating |2!j we focus 
here on methods that operate once the lattice potential 
has been raised. For example, one may sympathetically 
cool the atoms in the lattice using a second Bose-Einstein 
condensate 0, |3(J • A completely different approach is 
the filtering scheme of Rabl et al. • It operates in the 
no-tunnelling regime, transferring atoms between optical 
lattices so as to create a configuration with one atom 



per site. Such a loading scheme can, however, originate 
holes due to imperfections in the original cloud. Our 
analysis of filtering in the presence of a harmonic trap has 
shown that these holes, which are preferably located at 
the borders of the cloud, result in a considerable amount 
of entropy |32| ■ 

Here we propose and investigate several cooling 
schemes which overcome the limitations of filtering. The 
first set of schemes uses discrete operations to make 
atoms in different sites interact, thus concentrating the 
entropy on some atoms which are then expelled from 
the lattice. Due to the similarity with quantum infor- 
mation processing we term this kind of methods algo- 
rithmic cooling 46] . The second set of cooling methods 
combine filtering with either particle hopping or evapo- 
rative cooling techniques. All our cooling protocols are 
based on translation invariant operations (i.e. do not 
require single-site addressing) and consider the residual 
harmonic confinement present in current experiments. 
Although we will be mostly analyzing their effects on 
bosonic atoms, they can also be trivially generalized to 
fermions. 

This work has several spin-offs. First of all, we have 
designed algorithms that very efficiently remove all resid- 
ual defects from a Mott insulator, thus producing an en- 
semble of perfect registers for quantum computing 
Second, we propose how to create pointer atoms at the 
endpoints of a perfect quantum registers and show how 
these pointers can be used to tailor the register to a spe- 
cific length. Finally, since filtering is an important ingre- 
dient of all our algorithms, we have also designed a new 
incoherent version of filtering procedure which is better 
adapted to our protocols than the adiabatic scheme from 
Ref. [3ll. 

The paper is organized as follows. We start in Sect. iHl 
describing our system in terms of the Bose-Hubbard 
model and discussing the properties of current exper- 
imental states, such as densities, temperature and en- 
tropy. We also summarize the basic tools and figures 
of merit which underlie our cooling schemes. Sect, [ffi] 
is devoted to a particular tool, filtering. We describe 
the effect of filtering on the entropy and temperature for 
realistic initial states and study the optimal choice of ex- 
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perimental parameters. In the last part we propose an 
implementation of filtering which is similar to a radio- 
frequency knife. Once we have described all tools, in 
Sect. lIVI we introduce several protocols for cooling to the 
ground state. This includes the algorithmic cooling in 
Subsect. IIV Al and two other protocols in the following 
subsections. Based on similar ideas, Sect. |3 presents a 
quantum protocol that produces an ensemble of defect- 
free quantum registers and shows how to create pointer 
atoms on these registers. We conclude this work with 
some remarks concerning possible variants and combina- 
tions of our cooling schemes. In Appendix A we present 
a detailed description of the numerical method that we 
use to simulate classically correlated states. 



II. PHYSICAL SYSTEM 

A. Bose-Hubbard model 

We consider a gas of ultra-cold bosonic atoms which 
have been loaded into a three dimensional (3D) optical 
lattice. This lattice is created by six laser beams of wave 
vector k = 2tt/X propagating along three orthogonal di- 
rections. If the laser light is off-resonant with any atomic 
transition, the AC Stark effect induces a periodic poten- 
tial on the atoms of the form: 

V{x, y, z) = V 0x sm 2 (2kx) + V 0y sm 2 (2ky) + V Qz sin 2 (2fcz), 

(1) 

with a strength or "lattice depth" Vq proportional to 
the dynamic atomic polarizability and the laser inten- 
sity. The Gaussian profile of the laser beams creates 
an additional harmonic confinement which can be com- 
pensated by additional magnetic or optical confining el- 
ements [num. 

In the following we will mostly be concerned with one- 
dimensional (ID) lattices. In other words we will as- 
sume that the lattice potential is so strong along two 
directions, Vq v , Vq z ^> Vq x that tunnelling is only allowed 
along the third one. We will also assume that the con- 
finement along all directions is still much stronger than 
the atomic interaction strength. Under these conditions, 
the atoms can be described using a single band Bose- 
Hubbard model (BHM) [3^ , which for a lattice of length 
L reads 



Hbh 



L/2 

E 

fe=-L/2 



-J(a\a k+1 + /i.e.) + yn fc (n fc 



1) + bk n k 



(2) 

The parameter J denotes the hopping matrix element 
between two adjacent sites, U is the on-site interaction 
energy between two atoms and the energy b accounts 
for the strength of the harmonic confinement. Second 
quantization operators at and create and annihilate, 

respectively, a particle on site k, and = a\a,k is the 
occupation number operator. Since the tunnelling rate 
decreases exponentially with the trapping strength, Vq, 
while the on-site interaction U remains almost constant 



|34j , we have adopted this last value as the unit of energy 
for our work. Even more important, this exponential 
dependence allows one to reach the strong interaction 
regime U ^ J. At U « 11.6 J the ground state changes 
from a superfluid (SF) to a Mott insulator (MI) [H|3|. 
In the SF regime particles are delocalized over all lattice 
sites. In ID without harmonic trap the SF ground state 
reads: 



L/2-1 
,fc=-L/2 



|o>, 



(3) 



where N is the particle number. In the MI regime par- 
ticles are localized at individual lattice sites. For integer 
filling factors v — N/ L the ground state is given by: 



L/2-1 



MI 



-L/2 



(4) 



The presence of a harmonic trap typically leads to a wed- 
ding cake structure of MI states with different v. For very 
shallow traps one obtains only a single MI phase, which 
has filling v — 1 and which is centered around the bot- 
tom of the trap. In the Fock basis with respect to lattice 
sites this MI state can be written as: 

|'0o}mI = \0l/2 ■ ■ ■ 0_jv/2-ll-JV/2 • • • ljV/2-l0jv/2 • • • Ol/2-1/- 

(5) 

This state will be the target ground state for our cooling 
schemes. 



B. Initial states 

Throughout this paper we will work with ID thermal 
states in the grand canonical ensemble, which are char- 
acterized by two parameters: the temperature kT =1/(3 
and the chemical potential fj,. We are particularly inter- 
ested in the no- tunnelling limit 01, J — > 0. In this limit 
the Hamiltonian becomes diagonal in the Fock basis of 
independent lattice sites: {\n_ L / 2 ■ ■ ■ uq . . . n L / 2 -i}} and 
the density matrix becomes a product of thermal states 
for each lattice site, 



L/2-1 



p e 



Pk- 



(6) 



fe=-L/2 



This simplifies calculations considerably and for instance 
the von-Neumann entropy can be written as the sum over 
single-site entropies 



S(p)=tr(plog 2 p)=J2 S (Pk)- 



(7) 



Let us now study thermal states of the form (jSJ in 
more detail. In Fig. ^ we have depicted graphically all 
relevant energy scales in the no-tunnelling limit. The 
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FIG. 1: Illustration of important energy scales for bosonic 
atoms in an optical lattice in the no-tunnelling regime: on- 
site interaction energy U, chemical potential and harmonic 
confinement Vho = bk 2 . The characteristic size of the atomic 
cloud is given by fc M = J fi/b. 




FIG. 2: (a) Density distribution of two thermal states with 
equal Renyi entropy S 2 /N = 0.82 (S/N = 1 in MI phase) 
and equal particle number N = 65. At hopping rate J/U = 
(solid) the temperature is given by kT/U = 0.32 and at 
J/U = 0.16 (Vo = 5E r ) (dashed) one obtains kT/U = 0.46. 
The harmonic confinement is fixed at U/b = 370 (Vo = 22E r 
in transverse direction) . (b) The separation into two fermionic 
phases becomes clearly visible in the density profile of a ther- 
mal state at low temperatures. Numerical parameters: J = 0, 
kT/U = 0.072, fi/U = 1, U/b = 800, (jV = 65, s = 0.5). 



chemical potential determines the size of the cloud via 
the relation /i = bk^, where denotes the site at which 
= 0.5. For fi « U singly occupied sites at the 
border of the cloud become energetically degenerate with 
doubly occupied sites in the center. 

The analysis of recent experiments in the Mott regime 
[ttL l25j reveals a substantial temperature of the order of 
the on-site interaction energy U. This result is consistent 
with our own numerical calculations and translates into 
an entropy per particle s := S/N w 1. The particle num- 
ber in a ID tube of a 3D lattice as in ^(| typically ranges 
between N = 10 and N = 130 particles. A represen- 
tative density distribution corresponding to such initial 
conditions (with N = 65) is plotted in Fig. In this 
example the inverse temperature is given by (3U = 3.1. 
Since our cooling protocols lead to even lower tempera- 
tures, we will from now on focus on the low temperature 
regime, (3U S> 1. Moreover, we will only consider states 
with at most two particles per site, which puts the con- 
straint fj, < 2U — 1/P on the chemical potential. Such 
a situation can either be achieved by choosing the har- 
monic trap shallow enough or by applying an appropriate 
filtering operation |3lj . 

Under the assumptions e@ u ^> 1 and [i — U/2 > b + 
1/(2/3) we have shown in [32| that the density distribution 



of the initial state iJBJ can be separated into regions that 
are completely characterized by fermionic distribution 



1 



1 



/3(bfc 2 - M ) 



(8) 



To be more precise, at the borders of atomic cloud, bk 2 3> 
fi — U/2 + 1/ (2/3), the mean occupation number is given 
by (nk) ~ n\{k) with n\(k) := f k (b, (3, fi), while in the 
center of the trap, bk 2 <C — U/2 — 1/(2/3), one has: 
(n k ) ~ l + nu(k) with n\\(k) :— /fc(6,/3,/Xn) and effective 
chemical potential /in := fi — U. Note also that the 
underlying MI phase in the center of the trap is well 
reproduced by the function ni(k) , which originally has 
been derived for the border region. As a consequence, the 
density distribution for the whole lattice can be put in the 
simple form: (rik) ~ n\(k) +nn(k), which corresponds to 
two fermionic phases I and II , sitting on top of each other 
[Fig.Eb]. In other words, the initial state of our system 
can effectively be described in terms of non-interacting 
fermions, which can occupy two different energy bands 
I and II, with dispersion relations £1 = bk 2 and £n = 
bk 2 + U, respectively [Fig. 0. 

So far this fermionic interpretation is just based on 
the properties of the density profile of thermal states in 
the MI regime. In order to put this fermionic picture 
on more general grounds one can fermionize the BHM 
directly. In [32] we have shown that the dynamics 
at finite J can effectively be described by the following 
fermionic Hamiltonian: 



H 



cUk+i 



V2c\d k+1 + 2d\d k+1 + H. c. 



■;t, 



+ ]T [bk 2 c\c k + (bk 2 + U)d\d k 

k 



(9) 



Here, the fermionic operators Ck and d k refer to energy 
band I and II, respectively. This effective description 
is self-consistent as long as the probability of finding a 
particle-hole pair is negligible, i.e. (c k c\d\d k ) ps 0. We 
have shown above that this fulfilled for thermal states at 
low temperatures and at negligible tunnelling. Later we 
will see to what extent the validity can be extended to 
finite tunnelling rates. 

In general, the state of the system will however exhibit 
both classical and quantum correlations. In this case we 
have to represent the state in terms of a Matrix-Product 
State (MPS) [3^,|3^|. This concept can also be used to 
compute ID thermal states numerically |33|. Based on 
this method one can, for instance, estimate how the tem- 
perature of a ID tube changes when passing from the MI 
to the SF regime. Assuming that the process is thermo- 
dynamically adiabatic one obtains the new temperature 
with the following procedure. We tune the parameters 
of a thermal state in the SF regime until the expectation 
values for the entropy and particle number match the 
corresponding values of the initial thermal state in the 
MI regime. This can be illustrated with the following 
example. Starting from a representative state in the no- 
tunnelling regime with s = 1 and N = 65 we have to tune 
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FIG. 3: Effective description of thermal states in the no- 
tunnelling limit in terms of independent fermions occupying 
two energy bands. The dispersion relations are ei = bk 2 and 
£n = bk 2 + U, where k denotes the lattice site and U is the in- 
teraction energy. Increasing the harmonic trap strength from 
b to b' increases the chemical potential to fi' so that the pop- 
ulation of the upper band becomes energetically favorable. In 
the bosonic picture this process corresponds to the formation 
of doubly occupied sites. 



the figures of merit. For instance, the final entropy is 
given by S(p e ff). Given that our system can somehow 
thermalize, these will be indeed the properties of our fi- 
nal state. Finally, it is important to point out, that other 
variables, like energy or temperature, are not very well 
suited as figures of merit, because they depend crucially 
on external system parameters like the trap strength. 

Note also that in the case of quantum registers it can 
be erroneous to assume that a finite value of the final 
entropy implies the existence of defects. For example, 
we will propose a protocol below that generates a state 
which his an incoherent superposition of perfect quantum 
registers with varying length and position. This state 
has some residual entropy but it is ideal for quantum 
information processing. 



the temperature to kT = 0.46Z7 = 2.9 J at J/U = 0.16 
(Vb = 5E r ) in order to leave s and N unchanged [Fig.Efc]- 
Hence, in the SF regime one faces a substantial temper- 
ature of the order of the width of the lowest Bloch band. 
Note, however, that this is only a lower bound to the true 
temperature, because our approach does not include any 
sort of heating processes induced by the adiabatic evolu- 
tion. 



C. Entropy as figure of merit 

We will show below that algorithmic protocols are 
suited both for ground state cooling and the initialization 
of quantum registers. The goal in both cases is to create a 
pure state under the constraint of keeping a large number 
of particles. Given the fact that our protocols converge 
to the desired family of states, we can measure the per- 
formance of the protocol by computing the mixedness of 
the state. This mixedness can in turn be quantified using 
the von-Neumann entropy S Q. In some cases, such as 
for finite hopping J we will not be able to compute the 
von-Neumann entropy efficiently. We will then refer to 
the Renyi entropy 

S 2 = -log 2 (trp 2 ), (10) 

which is a lower bound S% < S and can be evaluated 
using MPS [Appendix A]. 

In order to assess the efficiency of a protocol in achiev- 
ing our objectives we define two figures of merit. The ra- 
tio of the entropies per particle after and before invoking 
the protocol, s//sj, quantifies the amount of cooling. The 
ratio of the final and initial number of particles, Nf/Ni, 
quantifies the particle loss induced by the protocol. Note 
that these figures of merit can sometimes be misleading 
and should therefore be applied with care. In the case of 
ground state cooling the entropy is only a good figure of 
merit if the state of the system after the cooling protocol, 
pt, is close to thermal equilibrium. If this is not fulfilled, 
we use an effective thermal state, pf — > p e ff, with the 
same number of particles, N, and energy, E, to compute 



D. Basic operations 

All our cooling protocols rely on a set of translationally 
invariant quantum operations that can be realized in cur- 
rent experiments with optical lattices. These operations 
are: 

(i) Particle transfer: Depending on the occupation 
numbers, an integer number of particles is transferred 
between internal states \a) and \b). This process can be 
described by the unitary operation 

U%+*' n - x :\m,n)~\m + x,n-x), (11) 

where x is an integer and \m,n) are Fock states with 
m and n atoms in internal states |a) and \b), respec- 
tively. Note that for this unitary operator it holds 
U™% x ' n ~ x = U™+ x , n -x- Certain operations, like C/ 2 °; 2 
or Ui'q , have already been demonstrated experimentally 
in an entanglement interferometer |38|. 

(ii) Lattice shifts: We denote by S x the operations 
which shift the \b) lattice x steps to the right. For 
example, S—i transforms the state ®jt|77ifc,njt}fc m to 
<8fc|rrc.fc, n.fc_|_i)fc. This operation can be realized in state- 
dependent lattices by adjusting the intensity and polar- 
ization of the laser beamsa 0, 0, HH 0, HH ■ 

(in) Merging and splitting of lattice sites: Making use 
of superlattices |4jj one can either merge adjacent lattice 
wells or split a single site into a pair of two sites. 

(iv) Emptying sites: All atoms in internal state \b) are 
removed from the system. We denote this operation by 
Ei,. It transforms the state ®k\mk,nk)k into ®k\mk,0)k- 
Experimentally, this can be achieved either by switching 
off the lattice potential acting on \b) or by coupling this 
state resonantly to an untrapped state. 

(v) Filtering: This means particle transfer operations 
of the form U m '™~ , followed by Eb, for all m > M. 
After tracing out the subsystem \b), the filter operation 
is described by a completely positive map acting solely 
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FIG. 4: Illustration of filter operation Fi which aims at pro- 
ducing a uniform filling of one atom per site. Defects arising 
from holes cannot be corrected. 



on atoms in state \a): 

F M [p\: ^p„, m |n)(m|^ (12) 

n,m 

^2 Pn,m\n)(m\ + 

\M)(M\. 



n,m<M 



n>M 



The first proposal for a coherent implementation of the 
operation Fi [Fig.0] appeared in Ref. [31f . More recently, 
a scheme based on resonant control of interaction driven 
spin oscillations has been put forward In we 

have proposed an ultra-fast coherent implementation of 
Fi relying on optimal laser control. In this paper we will 
propose an incoherent realization of filtering, which has 
the special virtue that it can be applied continuously. 




FIG. 5: Spatial dependence of the mean occupation number 
(n) and entropy S before (left) and after (right) the appli- 
cation of the filter operation Fi. The final particle distribu- 
tion can be well described by Eq. (JHJ (dashed line). Numer- 
ical parameters for initial state: Ni — 65, Sj=l, U/b = 700, 
fiU = 4.5 and fi/U = 1. Figures of merit: Sf/si — 0.56 and 
Nf/N, = 0.80. 



an energy range of 2kT around the chemical potential fi 
is given by: 



III. FILTERING 

We now study the filtering operation F\ in Eq. 1|12|) 
in more detail. This operation is especially relevant for 
cooling because it produces a state close to the ground 
state of the MI regime. Thus, it can serve as a benchmark 
which has to be beaten by alternative cooling schemes. In 
this section we begin summarizing the analytical results 
about filtering from Ref. j^J and discuss the conditions 
to reach optimal cooling efficiency. We close this section 
with a possible experimental implementation of filtering 
which is based on an incoherent coupling of atoms to a 
continuum of untrapped states. 



A. Theory 

In Fig. [5] we have depicted the particle and entropy 
distributions before and after one filtering step, F\. One 
observes that a nearly perfect MI phase with filling factor 
v = 1 is created in the center of the trap. Defects in this 
phase are due to the presence of holes and concentrate at 
the borders of the trap. This behavior is reminiscent of 
fermions, for which excitations can only be created within 
an energy range of order kT around the Fermi level. This 
numerical observation can easily be understood with our 
previous analysis of the initial state. Filtering removes 
phase II, which is due to doubly occupied sites, and leaves 
the fermionic phase I unaffected [Fig. |SJi] . 

The fermionic picture allows us to to find a simple 
estimate for the final entropy. From the expression for 
the density of states, g(E) = 1/^/bE, one immediately 
obtains that the number of states (lattice sites) within 



A := 



(13) 



This parameter is characteristic for the tail width of the 
density distribution [Fig. |S] and will in the following be 
central for the analysis of our protocols. Since the final 
entropy must be localized at these sites we expect that 
Sf ~ A. Indeed, the rigorous derivation in [32l | yields 
Sf = o\ A = o\ Nf/Pfi, with 01 = 7r 2 /(61n2). Here, we 
have introduced the final number of particles Nf, which 
can be well approximated by the parameter N: 



N 



(14) 



Hence, the final entropy per particle can be written as: 

1 



s 7 



Nf 



(15) 



which reflects the size of the region of defects A in units 
of the system size N. For the special choice /i = U (or 
equivalently (no) = 1.5) one finds the following expres- 
sions for our figures of merit |32| : 



Si 

Ni 




(16) 
(17) 



with numerical parameters on ~ 2.935 and rju — (1 — 
\ / 2)V^C(l/2) w 1.063. This result shows that filter- 
ing becomes more efficient with decreasing temperature, 
since s s /s l cx 1/VJ3U -> and N f /Ni -> 1 for (5U -> oo. 
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It is important to note that the state after filter- 
ing is not an equilibrium state, because it is energet- 
ically favorable that particles tunnel from the borders 
to the center of the trap, thereby forming doubly occu- 
pied sites. However, the system can easily be brought 
to equilibrium by adapting the trap strength. While 
tunnelling is still suppressed, one has to decrease the 
strength of the harmonic confinement to a new value 
with b' < bU / (2/i). The system is then in the equilibrium 
configuration /&(&', (3', //) © with rescaled parameters 
T = T b'/b and pi = p b'/b < U/2. This observation 
shows that it is misleading to infer the cooling efficiency 
solely from the ratio T'/T, because it depends crucially 
on the choice of b' . 
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B. Optimal initial conditions 

Let us now study how the cooling efficiency of filter- 
ing depends on the initial state variables 6, (3 and p. 
Since the initial temperature is dictated by the experi- 
mental setup, we consider only the trap strength b and 
the chemical potential p (which can be varied via the 
particle number ./V) as free parameters. Since our figures 
of merit are computed per particle we expect a weak de- 
pendence on iV and therefore focus on the b dependence. 
This can be studied in terms of the mean central occu- 
pation number (h)q. For instance, in the special case 
(no) = 1.5, we have obtained the analytical expression 
Ijlfijl for the ratio s//sj, which exhibits a strong temper- 
ature dependence. In contrast, for (no) = 2 one finds 
that the cooling efficiency becomes independent of the 
temperature: Sf/s,i w 1/V3 for (3U 1. This can be 
understood from the presence of a v = 2 MI phase in the 
center of the trap, which does not contain entropy. In 
the opposite regime (no) — 1 the protocol has no cool- 
ing effect at all. For general (no), we have computed the 
figures of merit numerically exact [Fig.^j. One observes 
that for initial entropies Sj < 1 a central filling (no) ~ 1.5 
is always close to optimal. 

A special situation arises when we approach zero tem- 
perature. As shown in Fig. Eli, the quantity Sf/si 
changes periodically when decreasing the trap strength 
b/U. In this regime an additional MI phase with v = 2 
is present in the center of the trap and the splitting of 
the entropy between the lower and upper fermionic band 
depends dramatically on the value of the harmonic con- 
finement. In particular, one can find trap strengths (e.g. 
U/b w 2000 in Fig. 01) at which the final entropy ap- 
proaches zero and only a comparatively small number of 
particles is lost (Nf/Ni = 0.9). For shallower traps the 
v = 2 MI phase starts to collapse. In this regime it be- 
comes very difficult to find proper thermal states which 
match the initial conditions in terms of entropy and par- 
ticle number. Finally, the central occupation approaches 
one and the cooling protocol leaves the initial state un- 
changed. 

It is also interesting to study the efficiency of filtering 



FIG. 6: Cooling efficiency of the filtering operation F\ for 
fixed initial entropy Si and particle number Ni as a func- 
tion of the harmonic confinement. The change of the trap 
is indicated by either the mean central occupation number 
(no) or the inverse trap strength U/b. We plot the figures 
of merit Sf/si (solid), Nf/Ni (dashed) and the weighted 
quantity (sfNi)/(siNf) (dotted). The initial parameters are 
N, = 100 and (a) s, = 1.5, (b) a. = 1, (c) s, = 0.5 and (d) 
Si = 0.05. 



when acting on the full 3D lattice. In this case only a 
small subset of ID tubes will satisfy the optimal initial 
conditions. Using the parameters of Sect. UTI fs,- = 1 and 
N = 2 10 5 ) we obtain that s f /s t = 0.78 and N f /N l = 
0.62. The relatively large particle loss results from the 
high densities in the center of the trap. This is also the 
main reason why the protocol performs worse compared 
to a ID tube with s, = 1 as shown in Fig. Eb- 



C. Experimental realization of continuous filtering 

Filtering constitutes a fundamental tool in all our cool- 
ing protocols. As shown before 0, 113, El filter opera- 
tions can be realized by coherently transferring particles 
between two internal states and then removing these par- 
ticles. What we will show below is that both processes 
can be combined producing an incoherent evolution that 
gives rise to the completely positive map F\ (I12JI . 




FIG. 7: Possible physical realization of incoherent and con- 
tinuous filtering. Atoms with energy E > U are resonantly 
coupled to a continuum of untrapped states and thereby re- 
moved from the system. 

The experimental procedure for achieving the map Fm 
is very simple. We will use two atomic states: one atomic 



7 



state shall be confined by an optical lattice, deep in a 
MI phase, while the other one will be in a continuum of 
untrapped states which are free to escape the lattice. We 
will couple the trapped and untrapped states with two 
Raman lasers which have a relative detuning of the order 
of the interaction energy, S ~ MU [Fig. Ej. As long as 
the coupling is active, the lasers will depopulate lattice 
sites which have too many atoms, n > M , while leaving 
other sites untouched |48j . 

If the untrapped states are such that they have few 
atoms and they are quickly expelled from the trap (for 
instance by a magnetic field gradient), these states will 
behave for practical purposes like a thermal bath in a 
vacuum state. If the coherence time of this bath, phys- 
ically determined by the time free atoms spend close to 
the trapped ones, is very short compared to the Rabi 
frequency Q of the Raman transition, we will be allowed 
to write a master equation for the trapped atoms. The 
solution of this equation converges at large times to the 
desired filtered state, e.g. M = 1 for F x . 

Conceptually, this mechanism is equivalent to the fre- 
quency knife from evaporative cooling, where atoms con- 
taining too much kinetic energy are expelled from the 
trap in order to lower the temperature. In our case, how- 
ever, it is the interaction energy we get rid off and, as 
a side effect, we make the filling of the lattice more uni- 
form. 

Compared to the optimal control scheme in 32], the 
operation of this much simpler method is not very fast. 
From the solution of the master equation it follows that 
states with occupation n > 1 decay after a time which is 
of the order of the inverse Rabi frequency, t% ~ l/(nfi). 
The main source of errors arises from the non-resonant 
coupling of the n = 1 state with the reservoir. The proba- 
bility of a defect (empty site) will approximately be given 
by: Po ~ ft 2 /U 2 . Hence, for po = 1CP 4 we get an oper- 
ation time t\ ~ 100/f which is comparable to the time 
scale of the adiabatic scheme jU- However, our inco- 
herent scheme has two big advantages. First, it can be 
applied continuously. Second, it does not put any con- 
straints on the interaction energies of the two species. 
This holds under the assumption that untrapped atoms 
are expelled so quickly from the trap, that they do not 
interact significantly with the trapped atoms. 



IV. GROUND STATE COOLING 

We have seen that the residual entropy after filtering is 
concentrated at the borders of the density distribution. 
Particles on these sites are also the only source of energy 
excitations because all doubly occupied sites have been 
removed. In the following we will propose several pro- 
tocols which selectively remove particles at the borders, 
thereby bringing the system closer to its ground state. 
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FIG. 8: Illustration of the algorithmic protocol. The state is 
initialized with the filter operation F%- (a) Operation U^'o- 
particles are transferred to state \b) (red), (b) Operation Sk ' 
Lattice |6) has been shifted ko sites to the right, (c) Iteration 
of the sequence: U±'®'q , E c , S-i. If a site is occupied by two 
particles of different species, then both particles are removed. 
Afterwards lattice \b) is shifted one site to the left. After k s 
iterations lattice \b) is shifted 2ko — k s sites to the left, (d) 
The same sequence as in (c) is applied, with the difference 
that lattice |6) is shifted to the right. The final distribution 
of atoms in state \a) is sharper compared to the initial distri- 
bution (dotted). Numerical parameters: Ni = 65, S2,j = 0.82 
{si = 1), U/b = 300, fc =42, k s = 20, N f = 30.2, s 2J = 0.19. 
After equilibration: S2j = 0.23 (s/ = 0.31). 



A. Algorithmic Cooling 

This is a cooling method based on a set of discrete 
steps which resemble a quantum information processing 
protocol. The steps are the following: (i) We begin with a 
cloud in thermal equilibrium in the no-tunneling regime, 
having two or less atoms per site, all in internal state |o). 
This can be ensured with a filtering operation F 2 . (ii) 
We next split the particle distribution into two, with an 
operation U 2 ' [Fig. [S^]- (iii) The |o) and \b) atoms are 
shifted away from each other until both ensembles barely 
overlap. We then begin moving the clouds one against 
each other, removing in the process atoms from doubly 
occupied sites. Experimentally, this can be achieved by 
introducing a third internal level |c) and applying the 
unitary operation U^'q in generalization of (|llfl . This 
sequence sharpens the density distribution of both clouds 
and it is iterated for a small number of steps, of order A 
[Fig. 0:]. (iv) The atoms of type \b) are moved again to 
the other side of the lattice and a process similar to (iii) 
is repeated [Fig. |8Ji]. (v) Remaining atoms in state |6) 
can now be removed. 

Notice that the final particle distribution will never be 
perfectly sharp [Fig. [SJi] because this protocol is intrin- 
sically limited by thermal and quantum fluctuations in 
the initial state. Qualitatively, in order to remove a par- 
ticle of type I a) from the tail, it has to coincide with a 
particle of type \b) from the other cloud. Errors arise 
when particles do not coincide and are thus proportional 
to the fluctuations of the density. If we want to clean 
about A sites or remove A particles, the errors will be 
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FIG. 9: (a) Density distribution after filtering (dotted) and 
after algorithmic cooling (solid) which can be well approxi- 
mated by a thermal distribution (circles) . For comparison we 
consider a variant of the algorithmic protocol, which exhibits 
no classical correlations at all times (dashed). Numerical pa- 
rameters: Ni = 100, s 3 ,i = 0.42 (si = 0.50), U/b = 1500, 
fc = 62, k s = 18, N f = 71.4, s 2 j = 0.0659; After equilibra- 
tion: s 2)f = 0.0740 (s f = 0.0988); Dashed line: k s = 28, 
Nf = 63.3, Sf — 0.0389. (b) Final entropy per particle 
Sf as a function of the initial particle number Ni for fixed 
Si = 0.5 (circles). For large Ni one obtains a 1/s/Wi de- 
pendence (solid). For comparison we plot Sf after filtering 
(diamonds) which is expected to be independent of Ni (dot- 
ted). Small variations in Sf can be attributed to the specific 
choice of the harmonic confinement. 



O(vA). In the limit f3U ^> 1, we obtain the following 
scaling behavior for the final entropy per particle 



A 



1 



1 



N 



(18) 



where N 114fl is the number of particles in \a) after step 
(ii) and Ni is the initial particle number. This simple 
estimate is already in good agreement with the exact 
result in |32(, which takes also particle losses into ac- 
count. Our findings show that the algorithmic proto- 
col becomes more efficient with increasing initial particle 
number. This is in contrast to filtering, for which s/ is 
independent of Ni . 

We have verified numerically the scaling (|T5|) [Fig. 13?]. 
The numerical simulation is by no means trivial because 
the protocol establishes classical correlations both among 
lattice sites and also among internal states. To repro- 
duce these correlations we have resorted to represent- 
ing the classical density matrices using MPS [Appendix 
A] . These numerical simulations show that the final den- 
sity distribution is very close to a thermal distribution 
[Fig. Indeed if we apply one or few iterations of 

the protocol, the final state will still contain some holes 
on the tails and be close to thermal equilibrium, as the 
computation of the Renyi entropy entropy S2j (|10|) be- 



fore and after equilibration shows [Fig.|5^]. On the other 
hand, for a large number of iterations, the final density 
matrix will be an incoherent superpositions of perfect 
uniform MI states which differ on their length and po- 
sition [Similar to Fig. U^J. These states are suitable for 
quantum computation but are far from equilibrium. 

It is clear that we could cool the atoms to the ground 
state with 100% efficiency if the density distributions of 
| a) and \b) atoms were perfectly correlated. It seems also 
that we can improve the performance in realistic situa- 
tions by removing the classical correlations which are es- 
tablished during the process [Fig. [3^,]. We also increase 
the performance if we break the correlations between the 
atomic species at the end of the protocol while leaving 
the inter-site correlations untouched. Each iteration of 
the protocol would then reduce the total entropy by a fac- 
tor of y/2 |32j| . From an experimental point of view, this 
means that the algorithm will work better when using 
multiple independent clouds to clean each other. These 
clouds may come from loading the lattice with atoms in 
different internal states, from splitting the lattice into 
multiple condensates, or simply by using the clouds from 
different ID tubes to clean each other. . . Many other pos- 
sibilities can be conceived. The protocol can be further 
improved by selecting specific particle number subspaces 
from the density matrix before restoring thermal equilib- 
rium (see Sect. 0for details). 

Note also that we need not get rid of the atomic cloud 
in state \b). For appropriate initial trap parameters, it 
is preferably to move this cloud back to the center and 
to pump all atoms back into internal state \a). The re- 
sulting state will exhibit a MI shell structure with two 
plateaus at densities n = 1 and n = 2. This is a simple 
way to engineer ground states of tighter traps. Finally, 
let us remark that this protocol will only correct inho- 
mogeneities along one direction. In a real experiment 
the protocol should be repeated along the other three 
directions. 



B. Non-algorithmic cooling schemes 

We have seen that algorithmic schemes perform much 
better than simple filtering. Here, we propose two non- 
algorithmic ground state cooling schemes which combine 
filtering with other quantum operations available in op- 
tical lattices. 



1. Sequential filtering 

The basic idea of this protocol is to first transfer parti- 
cles from the tails of the cloud to the center of the trap, to 
form doubly occupied sites, and then to remove these par- 
ticles with some filtering operation. To be more precise, 
one has to iterate the following sequence of operations: 
(i) Allow for some tunnelling while the trap is adjusted 
in order to reach a central occupation of (n)o ~ 1.5. 
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FIG. 10: Cooling efficiency S2,f/s2,i, in terms of the Renyi 
entropy, of filtering Fi as a function of the hopping rate J/U. 
For the numerical computation in terms of MPS we use a 
lattice of length L = 27, which typically contains N ~ 20 
particles. The trap strength is chosen to be b/U = 1/65 and 
/i = [7, yielding (n)o w 1.5. 



(ii) Suppress tunnelling and perform the filtering, (iii) 
Slightly open the trap so that final distribution resembles 
a thermal one (See Sect. 1111(1 . These steps are repeated 
sequentially until the entropy per particle converges. 

This protocol has been analyzed in detail in [3^1 an d 
the results can be summarized as follows. If one allows 
for equilibration then a state asymptotically close to the 
ground state can be reached after only a few filtering 
cycles. Without equilibration the final entropy is limited 
by the initial probability of finding a defect in the center 
of the trap. 

Alternatively, one can think of an implementation of 
sequential filtering that operates at fixed but non-zero 
hopping rate. Such a protocol would clearly profit from 
the fact that adiabatic changes of the lattice potential, 
which lead to additional heating, are not required. A key 
ingredient is the generation of doubly occupied sites in 
the center. One has to fix the hopping rate at a value, 
at which there exists a coupling between doubly occu- 
pied sites at the center and singly occupied sites at the 
borders. Our analysis of the BHM in terms of an ef- 
fective two-band fermionic model |32| shows, however, 
that this occurs typically only deep in the SF regime for 
J/U > 0.5. 

We have studied the sequential filtering at fixed J nu- 
merically, computing the cooling efficiency as a function 
of the hopping rate. The results are shown in Fig. ITU1 
For high temperatures the ratio sj/si is rather indepen- 
dent of J and we achieve some cooling. For very small 
temperatures, kT < 0.01/7, it changes dramatically with 
,/. This is a clear signature of the quantum phase tran- 
sition, which is expected to occur at J c w 0.085/7 in 
the thermodynamic limit. While for J < J c the state is 
well described in terms of independent wells and filtering 
works very efficiently, for J > J c particles are delocal- 
ized over the lattice, and filtering causes heating rather 
than cooling. Summing up, since we require a significant 
value of the hopping, J/U ~ 0.5, this second variant of 
sequential filtering can only be used as an initial step, 
when temperatures are still comparatively high. 



2. Filtering combined with frequency knife 

Particles located at the tails of the density distribution 
can also be removed with a method similar to evaporative 
cooling or a frequency knife. The idea is to make use of 
inhomogeneous on-site energies and to choose the detun- 
ing 6 of a radiation field in such a way that only atoms 
located at specific lattice sites are resonantly coupled to 
another internal state. Using a magnetic field gradient it 
has been demonstrated that individual lattice sites can 
be resolved within an uncertainty of about five sites |42| . 
In our case the spatial dependence of the on-site energies 
is naturally provided by the harmonic trapping potential. 
However, in order to make use of this inhomogeneity one 
has to ensure that the atoms are coupled to an internal 
state which responds differently to the ac-Stark shift in- 
duced by the lattice laser beams. Experimentally, this 
can, for instance, be achieved with a setup similar to the 
one for the creation of spin-dependent lattices Q, 0, |3!| . 
There, one atomic species is trapped exclusively by a~ 
polarized laser light, whereas the other species is trapped 
predominantly by a + polarized light. Hence, an optically 
untrapped internal state can simply be realized by using 
only <t + polarized laser light for creating the optical lat- 
tice. 

Let us now estimate the requirements on the the Rabi 
frequency of the transition depending on the trap 
strength b. For simplicity, we consider two internal 
states, and a configuration for which the excited state 
exhibits zero on-site energy at each lattice site. We are in- 
terested on removing particles which are typically located 
at a distance from the center of the trap. Hence, for 
resonant coupling of these atoms we must choose the de- 
tuning to be \S\ = bk 2 — \i. Particles at sites k^ + Ak feel 
an effective detuning S e ff(Ak) « bNAk, where N = 2k li 
denotes the particle number (|14|) . The probability that a 
particle at site fc M + Afc is transferred to the excited state 
is then given by: 

Pe(Afc) = n» + £W (19) 

In order to locally address the site k^ reasonably well one 
has to demand: f2 < bN. For typical harmonic trapping 
frequencies in the MI regime, uoto — y/8b/ (hmX 2 ) >2wx 
65 Hz and N = 50 this translates into fi < 1 kHz. 

Experimentally it should therefore be feasible to re- 
solve individual sites with an uncertainty of a few lattice 
sites and thus to sharpen the density profile within this 
uncertainty. Note that this scheme clearly profits from a 
large number of particles per tube. Another advantage 
is that it preserves the spherical symmetry of the density 
distribution in a 3D lattice. 

Finally let us remark that this method can easily be 
incorporated in the continuous filtering scheme [Fig. \7\ 
proposed in Sect. IIIII since there the atoms are also cou- 
pled to untrapped states. 
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ALGORITHMIC COOLING OF DEFECTS IN 
QUANTUM REGISTERS 



A perfect one-dimensional quantum register is a con- 
nected array of commensurately filled lattice sites. For 
most purposes the filling factor is v = 1. This state ap- 
pears naturally as the unique ground state of the BHM 
in the no-tunnelling limit and in the presence of har- 
monic confinement. Hence, any efficient ground state 
cooling protocol will produce a good quantum register. 
Moreover, we have seen in the previous section, that al- 
gorithmic protocols can be used to create an ensemble 
of quantum registers rather than a unique one (see also 
[M l32j|h In this section we will propose and analyze two 
alternative algorithms for the creation of a quantum reg- 
ister ensemble. As compared to previous proposals 0,113 
these schemes require only a small number of operations, 
which makes them more appealing for experimental im- 
plementation. The first protocol produces registers with 
filling v = 1, whereas the second protocol is optimized for 
filling v — 2. In addition, we propose how to transform 
atoms at the endpoints of each register into "pointer" 
atoms, which then enables addressing individual lattice 
sites. This can be used, for instance, to create registers 
of equal length. Even more important it offers the op- 
portunity to perform ensemble quantum computation in 
this system. 

We start from a rather cold cloud, (3U ^> 1, which 
has been subject to fundamental filtering operations Fm 
(|12[) . The result is an almost perfect MI in the center 
of the trap with some residual defects (or holes) which 
are predominantly localized at the borders [See Fig. [S] . 
Since we operate solely in the no-tunnelling regime these 
defects cannot redistribute nor evaporate. Our goal is 
to remove all these defects and we will achieve it by ap- 
plying nearest-neighbor quantum gates which simulate 
inelastic collisions between particles and holes. Simply 
put, whenever a particle sits next to a hole, the particle 
will be annihilated. This process is analogous to spin flips 
in ferromagnets and the formation of domains of equal 
magnetization. Thus, these algorithmic schemes can also 
be understood as controlled equilibration and cooling of 
defects. 

For all protocols that will be proposed below, the re- 
sulting state is a mixture of perfect (up to defects in the 
central MI phase) quantum registers, which differ only 
in their length and lateral position. The entropy of this 
state will be of the order Sf ~ log 2 A, where the parame- 
ter A <|13[) quantifies the translational uncertainty in the 
initial density distribution. Note, however, even though 
the final entropy is very small, these protocols are typ- 
ically not suited for ground state cooling. The reason 
is that the final state is far from thermal equilibrium. 
Numerically we find that the value of the entropy after 
equilibration is comparable to the value before invoking 
the protocol. However, if one makes use of the pointer 
atoms to select only registers with a specific length then 
these protocols indeed lead to cooling. 
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FIG. 11: Illustration of protocol 1 for the creation of quantum 
registers. After merging neighboring lattice sites, one empties 
single occupied sites. The algorithm stops, when all particles 
in the tails have been removed and only the central MI phase 



A. Protocol 1 

The operation sequence is illustrated in Fig. ^2 We 
first merge oddly aligned pairs of sites using a superlat- 
tice 0, El and empty sites with only one atom. The 
sites are again split and the original lattice is restored. 
This operation is then repeated several times, alternat- 
ing between even and oddly aligned sites, until the total 
entropy reaches a minimum. 

We have analyzed the performance of this protocol 
under realistic conditions using the MPS description to 
compute the at most classically correlated density ma- 
trix [Appendix A]. In Fig. El we plot the typical density 
distribution after different steps of the protocol. A single 
step changes very little the density but decreases dramat- 
ically the entropy (From S 2 = 17.5 to S 2 = 7.9). Indeed, 
the first steps account for the elimination of most defects. 
After a few iterations the value of the total entropy starts 
to saturate because the density matrix has collapsed to 
a classical ensemble of (almost) defect-free quantum reg- 
isters [Fig. IE]. As a consequence of the protocol, the 
number of atoms per register must be even, which leads 
to steps in the final density profile [Fig. I12j. 

Clearly the fix point of this protocol would be a state 
with no particles at all. Let us now estimate how many 
iterations M of the protocol are required to reach a 
state with reasonable particle number and tolerable de- 
fect probability. First, we have to point out that there 
are two sources of defects: (i) holes in the central MI 
phase and (ii) particles at the borders which are discon- 
nected from the central MI phase and which have not yet 
been erased by the protocol. The probability for defects 
of the first kind is negligible in the limit of low tempera- 
tures j3U ^> 1 [see Sect. Ill B| . Defects of kind (ii) can be 
assessed by the following observation: In order to erase a 
connected array of M particles, which is separated by at 
least one empty site from a central MI state, one requires 
exactly M/2 iterations of the protocol. The probability 
of finding defects after M iterations is then given by the 
probability of having states with array size larger than 
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FIG. 12: Left: Particle distribution after different steps of 
the protocol: after filtering Fi (dotted), after one iteration 
(dashed), after eight iterations (solid). Parameters of the 
initial state: N = 100, s = 0.7 (s 2 = 0.56), (7/6=1000, 
f3U = 5.8. Parameters after filtering: N = 70.5, s = 0.33 
(52 = 17.5, S2 — 0.25). Parameters after iteration 1: 
N = 65.5, S 2 = 7.89, s 2 = 0.12; iteration 2: N = 60.8, 
S 2 = 5.95, s 2 = 0.098; iteration 7: N = 48.8, & = 5.297, 
s 2 = 0.108; iteration 8: N = 46.8, S 2 = 5.296, s 2 = 0.113. 
Stopping the algorithm at iteration 3 gives a minimum in the 
entropy per particle: S2 = 0.095. Upon reaching thermal 
equilibrium this value has increased to S2 = 0.265 (s = 0.36) 
which is comparable to the value after filtering. However, 
when selecting the subspace containing only registers of length 
N — 48 then the protocol leads to cooling. The entropy per 
particle in thermal equilibrium will then be given by s = 0.17. 
Right: The step like structure of the density distribution al- 
lows one to deduce the states which contribute significantly 
to the final density matrix. The states can be classified ac- 
cording to their particle number and their lateral position. 
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FIG. 13: Illustration of protocol 2 for the creation of quan- 
tum registers. Spin-dependent lattice shifts allow to remove 
defects in a correlated way. 



peated this analysis for various initial conditions and our 
results confirm this observation. However, if one selects 
a particle number subspace with TV being larger than the 
mean value than our protocol can indeed be used for cool- 
ing. According to the data in Fig. ^| the entropy per 
particle can be reduced by roughly 50% as compared to 
filtering. 



M in the initial density matrix. This probability can eas- 
ily be computed numerically exact from the distribution 
(0. The derivation of closed expressions is however dif- 
ficult. Nevertheless, one can get a good estimate for the 
optimal number of iterations based on the following argu- 
ment. The characteristic tail width of the initial particle 
distribution JBJ is A H13|l . Hence, the occurrence of ar- 
rays of size M = 2A is already very unlikely. This implies 
that after roughly A iterations of protocol 1 we expect 
to have registers with negligible defect probability. This 
can be illustrated with an example. For the initial state 
in Fig. ^|we have A = 10 and the relative difference in 
the total entropy after the seventh and eighth iteration 
has reduced to 10 -4 , which implies a defect probability 
of the same order of magnitude. The typical size of the 
registers after A iterations is roughly N — 2A, where N 
(|14f> is the characteristic size of the initial Ml state. This 
implies N 3> A, otherwise no particles remain in the 
system. Since A N/PU, one has to require the low 
temperature regime (3U ^> 1. 

As mentioned above, the algorithmic ground state 
cooling protocol of the previous section can also be used 
for the creation of quantum registers. However, the 
ground state cooling algorithm involves O(N) operations 
as compared to 0(A) operations of the current protocol. 

Note also that the example presented in Fig. ^] indi- 
cates that this protocol does not lead to cooling. After 
restoring thermal equilibrium, the entropy has reached a 
value comparable to the one after filtering. We have re- 



B. Protocol 2 

We start from a state which contains only empty or 
doubly occupied sites. This can be achieved by applying 
filtering operation F%, followed by U^'q and Eb. The pro- 
tocol is depicted in Fig. ED We begin with transferring 
one particle per site to state \b). Then the | b) -lattice 
is shifted one site to the left and the same operation as 
before is performed. Afterwards one empties single oc- 
cupied sites. This procedure allows to remove defects in 
a correlated way. Occupied sites which have an unoccu- 
pied site to the right are emptied. Since the probability 
of finding particles in the center is close to one, the cen- 
tral MI is preserved except for losses at the right border. 
This procedure is repeated until all particle disconnected 
from the central MI phase become annihilated and only 
perfect MI phases in the center remain. One step of the 
protocol can be summarized in the following sequence 
of operations: U 2 ' Q , Si, U 2 'q, Eb, U^q, Eb- Following 
the discussion of protocol 1, this sequence has to be ap- 
plied approximately 2A times, where A denotes again 
the characteristic tail width of the initial particle distri- 
bution. The factor two stems from the fact that at each 
step of the protocol the size of "particle islands" in the 
tails, as well as the central MI phase, is reduced only by 
one as compared to two in protocol 1. The final den- 
sity matrix looks very similar to the one after protocol 1, 
with the difference that also registers with odd number 
of atoms appear. 
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C. Pointer atoms and register length control 

Given an ensemble of quantum registers we will now 
show how to create pointer atoms. To be more precise, 
our goal is to selectively transfer the two end atoms of 
each register to a different internal level. These pointer 
atoms enable single-site addressing which can used to 
engineer registers of specific length. 



1. Creation of pointer atoms 

Our scheme relies on the same set of operations that 
is used in current experiments for entangling atoms lo- 
cated at different lattice sites 0. We consider quantum 
registers with one atom per site. Initially all atoms are in 
internal state \a). A Hadamard transformation puts the 
atoms in the coherent superposition state (\a) + |6))/\/2. 
One then shifts the |6}— lattice one site to the right and 
waits an appropriate time until the on-site interaction 
between species |a) and \b) yields a collisional 7r-phase. 
This means that on each site the state \a)\b) is trans- 
formed into — |a)|6). After two lattice shifts to the left 
one waits again until a 7r-phase has built up. Then the 
lattice is shifted back to its original position and a second 
Hadamard operation is performed. The resulting state is 
a again a product state but the end atoms have been 
promoted to level \b). These atoms can be considered 
as pointer atoms because they mark the beginning (and 
the end) of each register and can thus be used to ac- 
cess every site within the register in a deterministic way. 
In practice, only one pointer atom is needed. The sec- 
ond pointer atom, e.g. the one on the left, can easily be 
removed by applying the following operation sequence: 
S-i, f/^i , Eb, U 2 ' ■ 



2. Manipulation of register length 

The pointer atoms can now be used to create an ensem- 
ble of registers of fixed length, a feature which is desired 
for quantum computation and which can even lead to 
cooling. We first present a protocol that requires only 
two internal states of the atoms. Then we show that the 
algorithm can be simplified considerably when a third 
internal level is included. In both cases we start from 
a situation where all registers have one pointer atom in 
state | b) which is located at the right most occupied site 
of the register. 

Protocol based on two internal states: The cen- 
tral idea is to remove atoms selectively from the system 
by transferring them to the pointer level \b). We first 
show how to discard all registers which are shorter than 

2 

a desired length M. We start with the sequence: U 1 \, 

S-ai, U ' 2 , Eb- This ensures that registers of length 
N > M are protected against further modifications, be- 
cause their pointer atoms have been removed. Next we 



promote atoms in doubly occupied sites to the pointer 
level and shift this (two-atom) pointer one site to the 
right. If the pointer hits an occupied site then one atom 
in | a) will be removed. Iteration of this process removes 
all atoms on the right of the current pointer position. 
To be more precise, one has to iterate M — 1 times the 

2 2 2 1 

sequence: U 2 'o, "1j ^0 2> ^i2> Eb- In a last step the re- 
maining doubly occupied sites are emptied via U^'q , Eb- 
After creating new pointer atoms the minimal register 
length is given by M' = M — 2. Again we keep only the 
pointer atom on the right side. Let us now show how 
to shorten registers of length N > M' to length M' . We 
first apply the sequence: S-i, U-y\, S-m>, S—\. For 
the target register of length M' this merely transfers one 
atom from the right end of the register to the left end. 
For larger registers one obtains a two-atom pointer that 
can now be used to select and discard all atoms located 
left of the pointer. This can be accomplished by iterat- 
ing the sequence: U^, Eb, U 2Q , S-± until all registers 
with appreciable weight in the density matrix have been 
shortened to the desired size. 

Protocol based on three internal states: This 
protocol is based on the operation Ui\' , which trans- 
fers an atom from \a) to \c), given that a pointer atom is 
present. This way one can use the pointer as a "marker" , 
which allows one to produce registers of desired length M 
in a very simple manner. One first marks an array of M 
atoms and then discards all atoms which have remained 
in level \a). The algorithm can be summarized as fol- 
lows: apply the sequence S-tM-i) ^U 1 \,Eb, U 2 q , U l \ ' ; 
iterate M — 1 times the sequence Si, ' ; finally apply 
the operation E a . 



VI. CONCLUSION 

We have proposed various ground state cooling 
schemes that allow to reduce the temperature in current 
optical lattice setups considerably. In particular, we have 
presented algorithmic cooling schemes which combine fil- 
tering with concepts inspired by ensemble quantum com- 
putation. We have shown that these algorithmic proto- 
cols can be designed both for ground state cooling and 
the creation of defect-free quantum registers. We have 
tested and analyzed our protocols under realistic experi- 
mental conditions, with special emphasis on the presence 
of a harmonic trap. We found that the residual entropy 
of all our protocols depends crucially on one parameter 
which accounts for the translational uncertainty of the 
initial particle distribution. 

A special virtue of our schemes is, that they rely on 
general concepts which can easily be adapted to differ- 
ent experimental situations. For instance, little modi- 
fications ensure that our protocols can be applied both 
to bosonic and fermionic systems. A second advantage 
of our protocols is, that they are designed for the no- 
tunnelling regime and hence do not necessarily require 
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equilibration processes induced by particle hopping and 
elastic collisions. Hence, they allow to approach the 
ground state of the no-tunnelling regime in very short 
times. 

In this sense, the collection of cooling schemes pre- 
sented in this article can be considered as a toolbox which 
is tailored for cooling atoms in optical lattice setups. The 
tool (or combination of tools) which is best suited for a 
given purpose, can be chosen according to the charac- 
teristic features of a specific experimental setup. For in- 
stance, in the case of large systems at high temperatures, 
one can think of combining filtering with the frequency 
knife method which is then followed by the algorithmic 
ground state cooling protocol. Or ground state cooling 
can be combined with adiabatic transformation of the 
Hamiltonian so as to produce ground states of models 
different from the simple Bose-Hubbard considered here. 
And one should also keep in mind that a 3D lattice struc- 
ture offers a large variety of possibilities, which have not 
been fully explored yet. 

The methods introduced here greatly enhance the po- 
tential of optical lattice setups for future applications and 
might pave the way to the experimental realization of 
quantum simulation and quantum computation in this 
system. We also hope that the concepts introduced in 
this work might trigger further research in the direction 
of ground state cooling in optical lattices. 
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APPENDIX A: NUMERICAL DESCRIPTION OF 
CLASSICAL DENSITY MATRICES IN TERMS 
OF MPS 

Our algorithmic protocols establish classical correla- 
tions between different lattice sites, when applied to ther- 
mal states in the no-tunnelling regime. Hence, a descrip- 
tion in terms of independent wells (as in JBJ) is no longer 
adequate. Therefore we refer to a representation in terms 
of matrix product states (MPS) |35ll3r|. 

To be more precise, for a ID lattice of length L we 
want to map a classical density matrix of the form 

P »' " 'I-..'/. ',...'/ (Al) 
W 

onto a pure state in MPS form: 

\iP) = Y,A\*AZ...A i £\i 1 ...i L ). (A2) 

w 

Here, A lk denote matrices of dimension D x D and 
ik G {0, . . . , d — 1} is the occupation number of site k. 



The matrices at the endpoints, A] 1 and A 1 ^ , are 1 x D 
and D x 1 vectors, respectively. The mapping from p 
to IV') can easily be accomplished by setting: p ll --- lL = 

Ail A12 A}L 

1 2 ' ■ ' L ' 

Expectation values for operators of the form O = 0\ <g> 
O2 ■ ■ ■ <E> Ol are calculated according to the relation 

L d 

(d)=trdp=l[J2^\ d kh)Al*. (A3) 

k=l i k = l 

Local operations on p, like unitary operations U™^ 1 
(|ll(l or filter operations l|12[l. amount to transformations 
of the local matrices AL* in the MPS picture. For illustra- 
tion, let us consider a completely positive map C which 
acts on a local state pk at site fc. In Kraus representation 
one can write 

C(p k ) = J2E aPk El (A4) 

a 

with the Kraus operators E a satisfying the completeness 
relation J2 a E l E a = 1- The MPS matr ices transform 
according to 

4 fc = EEi^ fe i 2 ^- ( A5 ) 

ot jk 

Non-local operations that involve more than one lattice 
site are more complicated to implement. As an example, 
let us consider the most complicated case, which occurs in 
our protocols: a spin-dependent lattice shift S—i which 
shifts the lattice \b) one site to the left. For this we 
have to consider two species of atoms. In generalization 
of (|X2|) . the MPS matrices A%' jk are now labelled with 
two physical indices, i k and jk, referring to states \a) 
and 1 6), respectively. It will turn out to be convenient 
to rewrite these matrices in tensor form: (A l £' 3h ) a ^ = 
A k (a,0,i k ,jk), with a,0 = X...D, 

To begin with let us consider only the first two lat- 
tice sites. We start out with tracing over subsystem \b) 
at the first site: Bi(a x , 0i, h) = A i («i, H, jl), 
because we want to omit sites which pass the system 
boundaries after shifting. The induced error is negli- 
gible, given that sites at the boundaries are not pop- 
ulated. We then multiply the matrix of the second 
site and perform the lattice shift: 9(ai, it, 02, ii, 32) = 
Bi(ai, 0,ii)A 2 (f3, f3 2 ,i 2 , j 2 ) -> Q{a x ,ii, 32,02,i2), where 
we use Einstein summation convention. After relabelling 
32 ji we perform a singular value decomposition: 

Q(ai,ii,ji, 02,12) = 

I7((ai,*i,ji),7) E(7,V) W^(j',(02,i2)), (A6) 

with E being a diagonal square matrix and 7,7' = 
1 . . . D < dD. The new MPS matrix for site 1 can then 
be defined as: 

(A\ lA ) :=U{ ai ,h,j 1 ,l)W 2 { 1 ,0 1 ). (A7) 
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Accordingly, we define 

B 2 {a 2 ,f3 2 ,i 2 ) ■■= S 1/2 (a 2 , 7 ) ^(7, i 2 ), (A8) 

and iterate this scheme until the end of the lattice. Note, 
that the dimension D of the MPS matrices can, in prin- 
ciple, increase exponentially with the number of lattice 
shifts. If D becomes larger than a desired value D max 
one has to resort to a truncation method similar to the 



one proposed for mixed quantum states |37| . However, in 
practice, we find that D increases only linearly with the 
number of lattice shifts, given that each shift is followed 
by a filter operation. Thus, with our method we can 
fairly easy simulate protocols exactly, i.e. without trun- 
cation, which involve up to 100 lattice shifts on lattices 
with up to L — 500 sites. 



D. Jaksch, H.-J. Briegel, J. I. Cirac, C. W. Gardiner, and 
P. Zoller, Phys. Rev. Lett. 82, 1975 (1999). [25 
J. Mompart, K. Eckert, W. Ertmer, G. Birkl, and M. 
Lewenstein, Phys. Rev. Lett. 90, 147901 (2003). [26 

G. K. Brennen, I. H. Deutsch, and C. J. Williams, Phys. 

Rev. A 65, 022313 (2002). [27 

0. Mandel, M. Greiner, A. Widera, T. Rom, T. W. [28 
Hansen, and I. Bloch, Nature 425, 937 (2003). 

K. G. H. Vollbrecht, E. Solano, and J. I. Cirac, Phys. [29 
Rev. Lett. 93, 220502 (2004). 
J. I. Cirac and P. Zoller, Science 301, 176 (2003). [30 
J. I. Cirac and P. Zoller, Phys. Today 57, 38 (2004). 
W. Hofstetter, J. I. Cirac, P. Zoller, E. Dernier, and M. [31 
Lukin, Phys. Rev. Lett. 89, 220407 (2002). 

J. J. Garcia-Ripoll and J. I. Cirac, New J. Phys. 5, 76 [32 
(2003). 

L.-M. Duan, E. Dernier, and M. D. Lukin, Phys. Rev. [33 
Lett. 91, 090402 (2003). 

J. J. Garcia-Ripoll, M. A. Martin-Delgado, and J. I. [34 
Cirac, Phys. Rev. Lett. 93, 250405 (2004). 

L. Santos, M. A. Baranov, J. I. Cirac, H.-U. Everts, [35 

H. Fehrmann, and M. Lewenstein, Phys. Rev. Lett. 93, 
030601 (2004). [36 

H. P. Biichler, M. Hermele, S. D. Huber, M. P. A. Fisher, 

and P. Zoller, Phys. Rev. Lett. 95, 040402 (2005). [37 

A. Micheli, G. K. Brennen, and P. Zoller, 
arXiv:quant-ph/0512222 [38 
R. Barnett, D. Petrov, M. Lukin, and E. Demler, 
arXiv:cond-mat/0601302 [39 
M. Greiner, O. Mandel, T. Esslinger, T. W. Hansch, and 

1. Bloch, Nature 415, 39 (2002). [40 

B. Paredes, A. Widera, V. Murg, O. Mandel, S. Foiling, 

I. Cirac, G. V. Shlyapnikov, T. W. Hansch, and I. Bloch, 
Nature 429, 277 (2004). [41 
H. Moritz, T. Stoferle, M. Kohl, and T. Esslinger, Phys. 

Rev. Lett. 91, 250402 (2003). [42 
T. Stoferle, H. Moritz, C. Schori, M. Kohl, and T. 
Esslinger, Phys. Rev. Lett. 92, 130403 (2004). 
M. Kohl, H. Moritz, T. Stoferle, K. Giinter, and T. [43 
Esslinger, Phys. Rev. Lett. 94, 080403 (2005). 

B. Laburthe Tolra, K. M. O'Hara, J. H. Huckans, W. D. [44 
Phillips, S. L. Rolston, and J. V. Porto, Phys. Rev. Lett. 

92, 190401 (2004). [45 
K. Xu, Y. Liu, J. R. Abo-Shaeer, T. Mukaiyama, J. K. 
Chin, D. E. Miller, W. Ketterle, K. M. Jones, and E. [46 
Tiesinga, Phys. Rev. A 72, 043604 (2005). 

C. Ryu, X. Du, E. Ye silada, A. M. Dudarev, S. Wan, Q. 
Niu, and D. Heinzen, arXiv:cond-mat/0508201 
G. Thalhammer, K. Winkler, F. Lang, S. Schmid, R. 
Grimm, and J. H. Denschlag, Phys. Rev. Lett. 96, 050402 



Uhrig, Phys. Rev. 
M. Troyer, 



and 



M. Kohl, and T. 



(2006). 

A. Reischl, K. P. Schmidt, and G. S. 
A 72, 063609 (2005). 
H. G. K atzgraber, A. Esposito, 

|arXiv:cond-mat/0510194l 

M. Kohl, |arXiv:cond -mat /0510567 
T. Stoferle, H. Moritz, K. Giinter, 
Esslinger, Phys. Rev. Lett. 96, 030401 (2006). 
T. Gericke, F. Gerbier, A. Widera, S. Foiling, O. Mandel, 
and I. Bloch, arXiv:cond-mat/0603590 
A. J. Daley, P. O. Fedichev, and P. Zoller, Phys. Rev. A 
69, 022306 (2004). 

P. Rabl, A. J. Daley, P. O. Fedichev, J. I. Cirac, and P. 

Zoller, Phys. Rev. Lett. 91, 110403 (2003). 

M. Popp, K. G. Vollbrecht, J. J. Garcia-Ripoll, and J. I. 

Cirac, Phys. Rev. A 74, 013622 (2006). 

F. Gerbier, A. Widera, S. Foiling, O. Mandel, T. Gericke, 

and I. Bloch, Phys. Rev. Lett. 95, 050404 (2005). 

D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and 

P. Zoller, Phys. Rev. Lett. 81, 3108 (1998). 

M. Fannes, B. Nachtergaele, and R. Werner, Comm. 

Math. Phys. 144, 443 (1992). 

S. Ostlund and S. Rommer, Phys. Rev. Lett. 75, 3537 
(1995). 

F. Verstraete, J. J. Garcia-Ripoll, and J. I. Cirac, Phys. 
Rev. Lett. 93, 207204 (2004). 

A. Widera, O. Mandel, M. Greiner, S. Kreim, T. W. 
Hansch, and I. Bloch, Phys. Rev. Lett. 92, 160406 (2004). 

G. K. Brennen, C. M. Caves, P. S. Jessen, and I. H. 
Deutsch, Phys. Rev. Lett. 82, 1060 (1999). 

S. Peil, J. V. Porto, B. Laburthe Tolra, J. M. Obrecht, 

B. E. King, M. Subbotin, S. L. Rolston, and W. D. 
Phillips, Phys. Rev. A 67, 051603 (2003). 

F. Gerbier, A. Widera, S. Foiling, O. Mandel, T. Gericke, 
and I. Bloch, arXiv:cond-mat/0601151 
D. Schrader, I. Dotsenko, M. Khudaverdyan, Y. Mirosh- 
nychenko, A. Rauschenbeutel, and D. Meschede, Phys. 
Rev. Lett. 93, 150501 (2004). 

T. Calarco, E. A. Hinds, D. Jaksch, J. Schmiedmayer, 
J. I. Cirac, and P. Zoller, Phys. Rev. A 61, 022304 (2000). 
J. Sebby-Strabley, M. Anderlini, P. S. Jessen, and J. V. 
Porto, arXiv:cond-mat/0602103 

L. Viverit, C. Menotti, T. Calarco, and A. Smerzi, Phys. 
Rev. Lett. 93, 110401 (2004). 

Note that our concept of algorithmic cooling of atoms 
has to be clearly distinguished from algorithmic cooling 
of spins which is a novel technique that allows to cre- 
ate highly polarized ensembles of spins in the context of 
NMR experiments, see e.g. P. O. Boykin, T. Mor, V. Roy- 
chowdhury, F. Vatan, and R. Vrijen, Proc. Natl. Acad. 



15 



Sci. USA 99, 3388 (2002). 
] The condition for the no-tunnelling regime in the pres- 
ence of a harmonic trap is given by |J/(6iV)j 2 <C 1, i.e. 
the hopping matrix element is much smaller than the 
average energy spacing between single particle states lo- 
cated at the borders of the trap. If desired one can also 
demand \J/b\ 2 <C 1, which ensures that single particle 
eigenfunctions even at the bottom of the trap are local- 



ized on individual lattice wells [45|| . 
[48] Additionally, particles from the tails may also be coupled 
to the continuum and expelled. In order to avoid that this 
process leads to an unwanted increase of entropy in the 
tails, one can lower the harmonic confinement once the 
MI phase has been reached. 



